Chromosome Research

The spatial organization of chromosomes is crucial for gene regulation and genome stability. Metaphase chromosomes appear physically discrete from one another in conventional chromosome preparations. However, interchromosomal connections, or "linkers", have been observed among mitotic chromosomes, and may play important roles in coordinated chromosome movements and genome stability. We use micropipette-based isolation and manipulation to analyze interchromosome linkers between mammalian mitotic chromosomes. Pulling an isolated chromosome reveals that thin linkers connect them; CREST staining indicates their locations to be at centromeres. The linkers display linear elasticity and have a length-doubling force of approximately 300 pN, comparable to the length-doubling force of an entire metaphase chromosome. Enzymatic treatments demonstrate that the linkers are disrupted by DNase, but not by RNase, protease, or microtubule polymerization (spindle fiber) inhibitors, indicating that linker connectivity is based on DNA. Immunofluorescence experiments indicate that histones and topo II are on the linkers, with topo II organized into discrete foci spaced by approximately 0.4 Mbp. After prolonged metaphase arrest using spindle inhibitors isolated genomes do not display interchromosome linkers. Finally, experiments on intact cells show that interchromosome linkers containing CENP-B can be observed between mitotic chromosomes, indicating they are not an artifact of genome isolation.

Cyclic episodes of chromosome compaction and de-condensation are features of eukaryotic cell division that aid mitotic segregation and help prevent aneuploidy. While biophysical data on mitotic chromosome structure has been previously obtained, heterogeneity within samples can confound analyses and precludes direct like-for-like comparisons. To circumvent this, we employed advanced flow cytometry to purify specific metaphase chromosomes with biotinylated telomeres from stably engineered mouse cells. We show that ESC-derived metaphase chromosomes 3 and 19 display distinct properties but share a conserved force-dependent mechanical response. In contrast, chromosome equivalents isolated from NSCs and preB cells show markedly different force-dependent responses, reflecting progressive differentiation stages. Covalent crosslinking of ESC-derived chromosomes alters biomechanical properties to mimic equivalents from more differentiated cells. Collectively, these results highlight the need to isolate specific, homogeneous metaphase chromosome samples to accurately decipher their complex behaviours.

Chromosome segregation is vital. Its disruption can cause aneuploidy, a hallmark of cancer. Centromere protein A (CENP-A) is an important protein during the segregation as it marks the location of the centromere, where the kinetochore assembles for microtubule attachment. Each mammalian centromere contains hundreds of CENP-A nucleosomes, whose spatial arrangement is expected to be critical for kinetochore function and error-free chromosome segregation. However, previous studies--mostly in fixed cells and focused on metaphase--have yielded conflicting results on CENP-A organization. Given the centromere's sub-diffraction size, we used stimulated emission depletion (STED) super-resolution microscopy to visualize the CENP-A organization in both living (about 45 nm resolution) and fixed (about 30 nm resolution) immortalized human retinal pigment epithelial cells (hTERT RPE-1). We found that CENP-A organization does not adopt a single architecture but spans a spectrum from dense clusters to fragmented subclusters, with mitotic stage-dependent abundance and morphology. CENP-A chromatin is most dispersed in prophase, and progressively compacts during prometaphase and metaphase as microtubules attach. In meta- and anaphase predominantly dense organizations are formed, often with a plate-like morphology. Yet, non-dense and non-plate-like organizations persist through metaphase and anaphase, suggesting CENP-A spatial reorganization is heterogeneous during cell division.

During meiosis, accurate chromosome segregation requires significant condensation and compaction. These processes are mediated by condensins, cohesins, and histone tail modifications. We identified that MET-2, a histone methyltransferase that catalyzes the dimethylation of histone H3 lysine 9 (H3K9me2), differentially impacts chromosome size in the male vs. female C. elegans germline. In met-2 null worms, autosomes during spermatogenesis are significantly larger than wild type, while chromosome size during oogenesis is unaffected. X-univalent size in males is also unaffected by loss of MET-2, indicating MET-2 differentially regulates autosomal and X-chromosome compaction in male spermatogenesis. Autosome size is not changed when males harbor a catalytically deficient MET-2 (met-2CD) or have mutations preventing germline histone H3K9 methylation (H3K9R). In addition, met-2 males, in contrast to met-2CD or H3K9R males, have more active RNA pol II in later stages of meiosis. These data suggest MET-2 plays a noncatalytic role in mediating chromosome structure and transcription. In met-2 male germ lines, genes on the X chromosome, which is typically enriched in H3K9me2, are significantly more likely to be upregulated than genes on autosomes, even though X-univalent size is unchanged. These results suggest that MET-2 plays a sex-specific role that is not limited to its enzymatic activity.

Knowledge of the nuclear organization of specific chromatin domains at the ultrastructural level remains limited. To address this, we employed correlative light and electron microscopy (CLEM) approaches to investigate the inactive X chromosome (Xi) in female human RPE1cells. Tandem fusion of macroH2A (mH2A), a histone variant enriched on the Xi, to eGFP and APEX2 enabled direct labeling and identification of the Xi by light microscopy, transmission electron microscopy, and electron tomography using conventional CLEM. The eGFP tag further allowed us to visualize the Xi at single nucleosome resolution under close-to-native conditions using cryo-CLEM. We found that, although heterochromatin domains of the Xi are substantially larger than those of autosomal regions, their chromatin density was not significantly different. Furthermore, our data reveal the organization of the interchromatin compartment (IC) within the Xi, characterized by a network of larger lacunae and smaller channels. Overall, we provide an ultrastructural characterization of a specific heterochromatin domain at unprecedented resolution, advancing our understanding of 3D chromatin organization and nuclear architecture in situ.

Sticklebacks (Gasterosteidae) encompass model organisms which are of particular interest for evolutionary and ecological genomics. Within Gasterosteidae, chromosome number is variable (2n=40-46) and independent fusions of homologous chromosomes have been proposed. The sea stickleback (or fifteen-spined stickleback, Spinachia spinachia ) has the lowest known number of chromosomes (2n=40) and hence is crucial in understanding chromosome evolution among sticklebacks, but is so far missing in genomic datasets. Here, we present a high-quality diploid genome assembly of S. spinachia. PacBio HiFi and Hi-C reads were assembled into a genome of 407.5 Mb in size, consisting of 20 chromosomes, with an N50 of 6.6 Mb and 98.96% complete single-copy BUSCO genes. A phylogenetic tree inferred across five stickleback species and four outgroup genomes from 19,156 genes, alongside synteny analyses and ancestral chromosome reconstructions, confirms S. spinachia as the sister species to the four-spined stickleback (Apeltes quadracus) and not as the sister group to all other sticklebacks as once thought. It has one species-specific chromosome fusion and shares two fusions with the three-spined stickleback (Gasterosteus aculeatus), none of which are present in its sister species. One of these fusions is also present in Pungitius, leading to reinterpretion of this fusion as ancestral to Gasterosteidae, with subsequent fission in Apeltes. This implies a lower ancestral chromosome number in Gasterosteidae (2n=44) than previously thought. The other fusion shared with G. aculeatus presents a case of convergence. Our results suggest that karyotype evolution in Gasterosteidae has been shaped by ancestral chromosome fusion, convergent fusion, and secondary fission.

Recent technological advances allow for the large-scale acquisition of genetic and morphological data: high-throughput sequencing has transformed the field of genomics while synchrotron X-ray microtomography enables rapid, noninvasive 3D imaging. However, integrating these approaches for the same specimens is challenging because X-rays can fragment DNA, and DNA extraction damages internal morphology, particularly relevant for small bodied organisms, such as insects. We systematically tested multiple extraction protocols and irradiation conditions across three model insect species. We irradiated more than 1,000 specimens under varying conditions and tested DNA quality through DNA barcoding and UCE sequencing. Our results demonstrate that high-quality DNA and high-resolution tomograms can be obtained from the same individuals, provided that the parameters are carefully optimized and rapid SR-CT scanning precedes DNA extraction. In this respect, our findings establish practical guidelines for combining genomics and phenomics, paving the way for comprehensive integrative digitization of biodiversity.

Poly(ADP-ribose) polymerases (PARPs) catalyze ADP-ribosylation, a conserved post- translational modification involved in DNA repair, transcriptional regulation, and chromatin remodeling. Although extensively studied in animals, the evolution and diversification of PARPs across the fungal kingdom remain largely unexplored. Here, we present the first kingdom-wide comparative genomic analysis of PARP proteins across 534 fungal species spanning eight phyla. We identified two primary conserved fungal PARP protein types corresponding to human PARP1 and PARP6. Both exhibit highly dynamic evolutionary histories characterized by frequent independent gain and loss events. Ancestral state reconstruction supports the presence of PARP1 in the last common ancestor of fungi, whereas the PARP6-like family has undergone repeated lineage-specific gains and losses. Fungal PARP6-like proteins retain a compact PARP catalytic domain fused to a C-terminal E2 ubiquitin-conjugating domain, whereas the PARP1 family displays extensive structural diversification through domain shuffling and lineage-specific fusions associated with DNA metabolism, chromatin remodeling, and signal transduction. Reconstruction of ancestral catalytic motifs across fungi and other eukaryotes revealed convergent evolution of a non-canonical H-Y-Y catalytic triad, with multiple motif variants co-occurring within individual proteins, suggesting functional diversification. In the Fusarium oxysporum species complex, we identified a lineage-specific expansion of the PARP family, driven exclusively by accessory chromosomes. Genomes with expanded PARP1 repertoires exhibited elevated basal PARylation, increased resistance to DNA-damaging agents that induce single strand breaks, and DNA damage- induced expression of accessory Parp genes. These findings reveal fungal PARPs as evolutionarily dynamic proteins that likely contribute to genome stability, adaptation, and pathogenicity.

Whole genome duplications leave lasting traces in our genomes. How these present in terms of gene content and order varies over time. While collinear blocks of paralogs, long stretches of conserved gene order and content termed microsynteny, are a distinctive feature of comparatively recent WGD and have been integral in reconstructing the history of ancestral duplication events, this signal degrades over time, making analysis of older events non-trivial. While gene order degrades quickly, gene content is often better conserved and recent work takes advantage of this to reconstruct older events and ancestral pre-WGD and post-WGD chromosomes. However, these new methods are complicated and not well-documented. Here we develop an automated and user-friendly pipeline for reconstructing ancestral chromosomes before and after WGD, and use the conservation of gene content to infer chromosomal rearrangement events in this timeframe. We verify the efficacy of our tool by reconstructing the ancestral acipenseriform, a model system for vertebrate WGD and rediploidisation. Our pipeline should serve to make ancestral reconstruction more accessible and provide a solid foundation for future analysis.

The cultivated dahlia has a complex evolutionary history that is likely the driver for its diverse floral forms. Previous studies have proposed cultivated dahlia as an allotetraploid, autotetraploid, and in one study an auto-allooctoploid. Here, we used PacBio HiFi long-reads and Dovetail Omni-C reads to assemble and annotate a genome for a tetraploid Dahlia variabilis cultivar Edna C. We identified Edna C as a segmental tetraploid with signature of auto- and allo-polyploidy events based on subgenome-specific repeat partitioning and single-copy gene tree construction within homoeologous groups. Additionally, we identified a cultivated dahlia-specific centromeric monomer shared between Dahlia variabilis Edna C and the Dahlia pinnata Kelvin Floodlight. Based on finding three major chromosome clusters, we propose that Dahlia variabilis Edna C is a segmental tetraploid with evidence of auto- and allopolyploidy.

Whole-genome duplication (WGD) causes chromosome instability through multipolar chromosome segregation driven by supernumerary centrosomes. WGD cells formed through distinct processes, mitotic slippage (MS) and cytokinesis failure (CF), show a prominent difference in viability after multipolar chromosome segregation: MS causes a more skewed homologous chromosome distribution than CF, resulting in more frequent nullisomic chromosome segregation with poorer survival through the first mitosis. However, the determinants of route-dependent differences in post-WGD cell viability remain largely unknown, particularly regarding the contribution of spatial rearrangement of supernumerary centrosomes. Here, we found marked differences in supernumerary centrosome distribution upon entry into the first mitosis after MS and CF, stemming from distinct nuclear geometry. The distinct centrosome distributions differentiated kinetochore capture patterning after MS and CF, whereas their modulations had minimal effect on the fidelity of subsequent chromosome segregation. In contrast, artificially extending the centrosome-centromere capture range by depleting the microtubule depolymerizer MCAK drastically suppressed the MS-linked aggravation of nullisomic chromosome segregation through equalizing chromosome capture by each supernumerary centrosome. These results suggest that centrosome-centromere capture range, rather than the spatial arrangement of the centrosomes themselves, determines the fidelity of chromosome segregation after WGD. Our findings provide fundamental insights into atypical cell proliferation mechanisms after WGD.

In the germline, small RNA-mediated gene regulatory pathways are associated with phase separated condensates known as germ granules across species. Germ granules can be challenging to purify based on their biophysical properties, therefore we used an in vivo proximity labeling approach, TurboID, to define the proximal interaction networks of the germ granule-localized Argonaute proteins: CSR-1, PRG-1/Piwi, PPW-2/WAGO-3, WAGO-1 and WAGO-4, and the nuclear-enriched Argonaute HRDE-1 as a comparator, in C. elegans. We found that although CSR-1 and WAGO-4 bind to the same set of 22G- small RNAs and target an overlapping set of germline genes, they possess distinct proximal interaction networks. Displacement of WAGO-4 from the germ granules highlights distinct interactors involved in chromosome biology and points to a potential nuclear role for WAGO-4. PRG-1/Piwi, WAGO-1, and PPW-2/WAGO-3, act in a common pathway and share an interaction network, while HRDE-1 displays a distinct interaction network that highlights its potential to influence splicing, rather than its widely accepted role in histone modification. Performing TurboID on other sRNA pathway factors, including the RdRP RRF-1 led us to uncover a potential nuclear role for this protein, and overlapping our TurboID data with published proximity labeling data enabled us to identify a core set of eleven mostly uncharacterized proteins that play roles in fertility, germ granule organization, and transgenerational epigenetic inheritance. One of these factors, D2005.4, appears to be a heme-binding germ granule-localized protein with a Tetratricopeptide Repeat Domain that interacts with all germ granule AGOs and may play scaffolding or enzymatic roles in germ granules. Collectively, our work highlights the utility of proximity labelling to uncover insights into AGO and sRNA pathway factor function and identify new factors that could be of importance in germ cell biology.

Targeted protein turnover of the transcription factor C/EBP is enhanced by the Tribbles (Trib) pseudokinase family of adaptor proteins, which bind to a conserved degron motif in C/EBP substrate and recruit a cognate COP1 E3 ligase to ubiquitinate C/EBP for degradation. We identified a conserved degron in the Drosophila C/EBP ortholog Slbo and fused this motif to a fluorescent protein (GFP or RFP) under the control of a constitutive promoter. Transgenic flies expressing this fluorescent protein (FP)::degron transgene, which we call the Fluorescent Ubiquitination-based Trbl Turnover Indicator (or FUTTI) sensor, showed reduced levels of fluorescence in response to increased levels of Drosophila Tribbles (Trbl) and conversely, increased levels when Trbl activity was reduced. Mutations in the Trbl conserved pseudokinase substrate binding pocket and in a divergent PEST-like motif in the C-terminal tail region disrupted Trbl activity in the FUTTI assay. These data indicate that FUTTI will be a valuable tool to identify sites of Trbl activity in vivo and the critical features of the protein required for function.

Dinoflagellate chromosomes adopt a highly condensed and organized morphology, with periodic bands and arches observed by traditional Transmission Electron Microscopy (TEM). However, the limited two-dimensional field of view of TEM has prevented a precise characterization of the inherently three-dimensional organization of dinoflagellate chromosomes. Moreover, given the vast diversity among dinoflagellate species and the lack of a systematic comparison of their chromosomes, it remains unclear whether dinoflagellate chromosomes share common organizational features or instead exhibit significant cell- or species-specific differences. Here, we acquire three whole-nucleus 3D Focused Ion Beam Scanning Electron Microscopy (FIB-SEM) datasets at 4 nm voxel size for each of four dinoflagellate species: Symbiodinium microadriaticum, Breviolum minutum, Fugacium kawagutii, and Crypthecodinium cohnii. We compile these data with previously published image volumes from four additional species and present an analysis of the largest collection of dinoflagellate FIB-SEM images to date. Common features observed across all eight species include the absence of physical confinement or spatial clustering of chromosomes in the nucleus. In addition, by decomposing each chromosome into a weighted sum of orthogonal shapes using Spherical Harmonics Expansion, we find a principal component encapsulating 88% of the total shape variance that is common to all species. However, our analysis also reveals differences in chromosome morphology across species. First, while many chromosomes exhibit surface ridges with left-handed helical twist, the proportion of chromosomes with such ridges varies extensively across species. Second, while chromosomes in most species are discrete and well-separated, chromosomes in F. kawagutii are interconnected in a single contiguous network. Lastly, to our knowledge, we report the first observation in eukaryotic cells of toroid-shaped DNA objects, whose numbers vary dramatically across cells and species. Overall, our results show that dinoflagellate chromosomes exhibit both shared organizational features and pronounced species-specific deviations.

Insect pigmentation is an ecologically critical trait influencing many physiological processes. In Drosophila melanogaster, abdominal pigmentation is sexually dimorphic: males have fully pigmented posterior segments, while females exhibit a posterior melanin stripe. Pigmentation relies on the expression of pigmentation genes that encode enzymes involved in pigment synthesis. These genes are tightly regulated during pupal and young adult stages. To expand the gene regulatory network of pigmentation genes, we conducted an RNAi screen using the yellow-Gal4 driver, expressed during the pupal stage in abdominal epidermis. One of the candidates from this screen, Kdm5, encodes a histone demethylase erasing the H3K4me3 histone mark catalyzed by the histone methyl-transferase Trithorax (Trx). We show that Kdm5 down-regulation reduces abdominal pigmentation, mimicking trx down-regulation. Kdm5 activates melanin production through regulation of the pigmentation gene tan. Transcriptomic analyses reveal that Kdm5 and Trx share many targets in pupal abdominal epidermis, including piRNA pathway components such as piwi and aubergine. These piRNA components, originally associated with transposon silencing in the germline, also function in some somatic tissues such as the nervous system, the fat body or the gut. We demonstrate that Piwi and Aubergine participate in female abdominal pigmentation establishment, without evident piRNA production. We also show that Kdm5 and Piwi act not only in pupal abdominal epidermis but also in pupal fat body. This study therefore expands the regulatory network of pigmentation genes. It identifies a new somatic function for Kdm5 and Piwi and reveals a role for pupal fat body in female abdominal pigmentation regulation.

A complete reference genome for species with chromosomally-determined separate sexes should contain scaffolds for all sex chromosome homologs. However, sex chromosomes present distinct computational challenges compared to autosomes. Here we present a k-mer based analysis that utilizes whole-genome sequencing of a few sex-identified isolates: Cytogenetics-By-Sequencing (CBS) tools. Unlike other approaches that typically address one aspect of the sex chromosomes, CBS-tools strives to guide users from the discovery of the heterogametic sex through identifying the sex-determination region (SDR). The core of CBS-tools is automated quantification of sex-specific k-mers in order to predict the heterogametic sex. Using publicly-available datasets, CBS-tools correctly identified the known sex-system of the 31 species tested. Additionally, we used these k-mers to verify and correct phasing of sex chromosomes between haplotypes in species representing different sex-systems. Finally, we used these k-mers to delimit the SDR boundary using an interactive web platform. CBS-tools was developed with previously unexplored sex chromosome systems in mind, but is also suitable for well-examined sex chromosome pairs.

Satellite DNA is long arrays of tandem repetitive DNA located often near the centromeres of chromosomes, whose function, or lack of, has been debated since its discovery. Although situated in heterochromatin, satellite DNA may be expressed as long noncoding RNAs (lncRNAs). Although there are a few examples of satellite lncRNAs being characterized, and functions suggested, how widespread and functionally important they may be for developmental processes is not understood. Here, we take an evolutionary approach to investigate satellite lncRNA expression in Drosophila spp. ovaries, a tissue whose development is well-characterized but where satellite expression has only been minimally explored. Using a publicly-available total RNAseq dataset, we find that 118/156 surveyed satellite DNAs were expressed across 10 species, with 33 satellites having high expression over 20 RPM. However, all but two of these expressed satellites (AAACTAC in D. virilis and ACAGACAGACAGG in D. ananassae) had higher read counts in a sister smallRNA dataset, suggesting that most satellite transcripts primarily serve as precursors for piRNA biogenesis. The two "stand-alone" lncRNAs were highly strand-biased, with 96-97% of the total reads coming from one strand. We further investigated AAACTAC expression with RNA FISH and found the transcript is specifically present in the oocyte nucleus following a dynamic spatiotemporal pattern, with the highest expression in stage 3-5 oocytes. The transcription pattern of AAACTAC is conserved in the three other virilis clade species that contain this satellite DNA. Further, we found expression of unrelated satellites in more distantly related D. borealis and littoralis both in the oocyte and the nurse cells. Overall, our work identifies a novel lncRNA AAACUAC found in the early oocyte nucleus, which is conserved across ~5 MY of evolution, and is therefore a strong candidate for the discovery of novel functions of satellite lncRNAs in development.

Bipolar spindle assembly is essential for accurate chromosome segregation. KIFC1, a conserved Ran- regulated minus-end-directed kinesin-14 motor, accumulates in the nucleus during interphase and promotes chromatin-mediated spindle assembly during mitosis and meiosis. In human oocytes, reduced KIFC1 levels destabilize meiotic spindles, a defect that can be rescued by increasing KIFC1 expression. However, how KIFC1 expression levels affect mitotic spindle stability during cleavage divisions in vertebrates remains unclear. Here, we show that whereas an approximately 50% reduction in KIFC1 causes no detectable defects in spindle assembly, approximately 10-fold overexpression of KIFC1 induces monopolar spindle formation, leading to chromosome mis-segregation and embryonic lethality in medaka early embryos. KIFC1 overexpression results in ectopic centrosomal localization during interphase, impairing the separation of duplicated centrosomes before mitotic entry. Analyses of KIFC1 mutants demonstrated that these centrosome separation defects require KIFC1s microtubule-binding and motor activities and are further enhanced by deletion of KIFC1s nuclear localization sequences. Together, our findings demonstrate that tight regulation of KIFC1 expression and its nuclear sequestration is essential for the proper separation and positioning of duplicated centrosomes before mitotic entry, thereby ensuring efficient bipolar spindle assembly during the rapid cleavage divisions of vertebrate embryos. HighlightsO_LIKIFC1 accumulates in the nucleus and at the embryonic spindle midplane via the Ran pathway. C_LIO_LIPartial KIFC1 depletion does not impair spindle assembly in medaka early embryos. C_LIO_LIKIFC1 overexpression induces monopolar spindles by preventing centrosome separation. C_LIO_LICentrosome separation defects require KIFC1 microtubule-binding and motor activity. C_LI

Polyploidy can be a critical factor for explaining plant trait variation, niche diversification, or speciation. However, inferring ploidy from silica-dried or historical samples using chromosome counts or flow cytometry is not possible, and scaling up ploidy estimation to population-level fresh contemporary samples can be challenging as well. Thus, we present a new method for estimating ploidy levels directly from sequencing data using machine learning; the Polyploid Population Genomics Tool Kit (PPGTK). The machine-learning approach is advantageous as it relaxes the assumptions of previous probabilistic methods and provides per-sample probabilities, allowing investigators to evaluate uncertainty in their system of interest.. We demonstrate performance and accuracy of the method on simulated and empirical data. Simulations showed above 99% accuracy, even for low coverage data, as long reads were mappable to the reference genome. For empirical analyses, we used target enrichment data from blueberry wild relatives (Vaccinium sect. Cyanococcus) and whole-genome data from sweetpotato wild relatives (Ipomoea ser. Batatas). Ploidy was recovered with 99% accuracy across 70 Vaccinium individuals and 97% across 82 Ipomoea individuals. Analysis of many individuals is fast and requires only a multisample VCF, which is presumably generated for the research anyway, and some samples of known ploidy for training the classifier. The approach implemented in PPGTK is promising for collections-based research as well, enabling ploidy classification of historical specimens based on present-day observations. The method is implemented in a new Python package as a single command that can run on a conventional laptop.

Primordial germ cells (PGCs) are the stem cell lineage responsible for transmitting genetic information across generations and are therefore central to studies of germline development, reproduction, and evolution. However, a lack of molecular tools in jawless vertebrates such as lampreys and hagfish has limited comparative insights into the evolutionary origins of vertebrate germline cells. Here, we report the development of the first germline-specific reporter (GSR) in the sea lamprey (Petromyzon marinus), Tg(piwil1: egfp-UTRnanos1), a piwil1-based germline reporter that is driven by the sea lamprey piwil1 regulatory region that directs Green Fluorescent Protein (GFP) expression in early PGCs. Promoter activity was validated by robust GFP expression in PGCs of zebrafish embryos co-injected with Tol2 transposase mRNA. Injection of the construct into sea lamprey embryos resulted in early, persistent, and gonad-localized GFP-positive cells. These results demonstrate stable germline-specific transgene expression in a jawless vertebrate and suggest deep evolutionary conservation of germline regulatory mechanisms between jawless and jawed vertebrates. Furthermore, the piwil1-based germline reporter line enabled lifetime visualization of zebrafish germline cells, revealing that increased PGC abundance during embryogenesis biases toward female sexual differentiation. Together, this study establishes the first GSR system in a jawless vertebrate and provides a useful platform for investigating germline biology and vertebrate germline evolution.